Nodal domains in open microwave systems 
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Nodal domains are studied both for real tpR and imaginary part of the wavefunctions of an 
open microwave cavity and found to show the same behavior as wavefunctions in closed billiards. In 
addition we investigate the variation of the number of nodal domains and the signed area correlation 
by changing the global phase tp g according to ipR + tyl = e t{Pg (ipR + i^'i ) • This variation can be qual- 
itatively, and the correlation quantitatively explained in terms of the phase rigidity characterising 
the openness of the billiard. 
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I. INTRODUCTION 

Approximately 200 years after Chladni's work on 
acoustic plates detecting and classifying nodal lines [l| 
there is a revival due to two pioneering works on the sub- 
ject of nodal domains. The first to mention is by Blum 
et al. [IJ proposing nodal domains as a tool to separate 
integrable and chaotic systems from each other. They 
explained their numerical findings for nodal domains in 
chaotic billiards in terms of the random plane wave ap- 
proach 0] and found a good agreement. Bogomolny and 
coworkers 0] introduced a percolation model, yielding 
explicit results, e. g., for the linear increase of number of 
nodal domains with the quantum number with a slope of 
0.0624. These results have been experimentally verified 
in closed microwave billiards 0, [|| . 

Other studies of nodal domains are concerned with the 
continuation into the classically forbidden region taken 
tunneling into account 0], with graphs Q, with iso- 
spectrality on graphs @ or billiards, i.e., whether the 
shape of the drum can be determined by nodal domain 
counting [n| [HI, IHj] > with quantum maps [HI , and with 
the anharmonic oscillator fl4| . 

Here we will concentrate on real physical systems, 
which are always open due to the measurement process 
and absorption. We shall show that nodal domains still 
can be defined and have the same behavior as in closed 
systems. In addition it will be illustrated that the vari- 
ation of the number of nodal domains in dependence of 
a global phase of the wavefunction can be used as a tool 
to determine the 'openness' of the billiard. 



II. NODAL DOMAINS FOR REAL AND 
IMAGINARY PARTS 

The billiard used is a quantum dot like structure with 
two attached leads of width 2 cm (for more details see 




FIG. 1: (color online) Real part tpn (left), imaginary part ipi 
(center), and modulus 2 = ipit+ipi of the wavefunction ip at 
v = 5.64 GHz corresponding to a Weyl number of nweyi=32.5. 
Correlations between if)R and tpi have already been removed 
by a global phase, see text. Additionally the nodal points of 
the modulus and the nodal lines of real and imaginary parts 
are marked. Nodal points occur at crossings of a nodal lines 
of %pE, with that of ipi. 

[HI]). Additionally we added two insets to two of the 
sides (see also Fig. [TJ to avoid any bouncing ball struc- 
tures jl5l . HH ] . The measurements have been performed 
on a grid of step size 2.5 mm. A similar billiard has al- 
ready been used to determine long range correlations in 
the wavefunctions [r^ | and vortex distributions and cor- 
relations [l8| . For each point the transmission 612 from a 
fixed antenna 1 in one lead to the probing antenna 2 has 
been measured including the phases. As the transmission 
depends only on the field distribution at the positions of 
the antennas, and the position of antenna 1 is fixed, one 
can obtain the electric field distribution E z (x,y) at the 
position of antenna 2 0,E3]- We will use the term 'wave- 
function' for this quantity as it is related to the quantum 
mechanical problem of a particle in a box, even though 
it is not the eigenfunction of a single resonance but a 
superposition of different cigcnfunctions. From this mea- 
surement the wavefunction ip(x,y) = E z {x,y) has been 
determined including its real and imaginary parts (see 
[H, HH and chapter 6 of (H). In Fig. [T] real part i/)r, 
imaginary part vjj and the modulus of a typical wave- 
function is shown. 

Fig. [^left) shows a plot of the imaginary part vj'j ver- 
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FIG. 2: Imaginary versus real part of the wave function 
at v = 13.84 GHz, for the directly measured wave function 
(left), and after uncorrelating them (right) by ipR + iipi — 
e -"p g ,o (^p' R -)-j^), ip g i s a global phase from the experi- 
ment acquired due to the antenna and leads. 



sus the real part ip' R for a typical wavefunction prior to 
any correction. Each dot corresponds to a measured grid 
point. Usually the cloud of dots is not aligned with its 
main axes with respect to the coordinate axes leading to 
a correlation of real and imaginary part. This correlation 
has its origin in unwanted global phase shifts ip 9t Q mainly 
from the lead and the antennas. This global phase has 
been removed together with the correlations between ipR 
and ipi by means of a phase rotation, 
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where > ("4>j)i thus adjusting the main axes of the 
cloud of dots to the coordinate axes (Fig. [2] (right)). For 
a completely open system a circular cloud would have 
been expected. The eccentricity thus reflects the lack 
in openness (for this wavefunction). For a quantitative 
description of the openness the phase rigidity \p\ 2 has 
been introduced [HJ . p is defined by 



(1% 



(2) 



This quantity is immediately available from the cloud 
of dots as shown in (Fig. fright)). It is noteworthy to 
mention that the phase rigidity is a highly fluctuating 
quantity and is different for each wavefunction. 

Because of the openness of the billiard there is a strong 
overlap of eigenfrequencies at each frequency making a 
direct determination of the mode number n impossible. 
Therefore the Weyl number [23| has been used to deter- 
mine this quantity, in our billiards given by 
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C 



(3) 



where fc = 2ttp/c is the wavenumber, v the frequency, c 
the speed of light, A and S are the area and the circum- 
ference of the billiard and C is a constant determined by 
corners and the curvature in the billiards. The order of 
the constant is 1 and in our case we used A as 0.0361 m 2 , 
S re 0.807m and C=0. Strictly speaking the Weyl for- 
mula is not applicable for open systems, since the res- 
onances are shifted from the real axis into the complex 
plane. Occasionally resonances even may be removed 
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FIG. 3: (color online) The number of nodal domains u n ver- 
sus the Weyl number nwcyi for the real (crosses) and imagi- 
nary part (diamonds) . The red dashed line correspond to the 
theoretical prediction v n = 0.0624n of the percolation model 
Q. The blue solid lines are fits including boundary effects @] 
with Eq. ((6} with a linear slope of 0.059 and 0.060 for real 
and imaginary part respectively. 



from the spectrum due to a strong coupling to the sur- 



roundings 



24J , and fractal Weyl laws may occur |25|, |26j . 
But still the Weyl formula is the best approximation 
known in this case to obtain the mode number n. 

From the measured real and imaginary parts we de- 
termined the nodal lines by means of a bilinear interpo- 
lation. From the bilinear interpolation we get the lines 
within the grid. Using this information and the Hoshen- 
Kopelman method [27J we obtain the number of nodal 
domains v n and their areas. 

In the case of open systems the distribution of inten- 
sities P(|V>| 2 ), e.g., can be obtained from the random 
superposition of plane waves (RSPW) approach Q, ex- 
pressing ip as a 



#0 = I>(k)< 



(4) 



In open systems the plane wave amplitudes <x(k) have a 
Gaussian distribution with zero mean and with variance 



(a(k)a(-k)) = p(o(k)o*(k)), 



(5) 



where \p\ 2 is the phase rigidity as given by Eq. (J2]). In this 
model the phase rigidity is only defined by an average and 
individual realizations therefore have fluctuating phase 
rigidities, depending on the number of plane waves used. 

First we investigate the number of nodal domains v n 
for the real and imaginary part. Bogomolny and Schmidt 
[H predicted from the percolation model a linear increase 
with a slope of a=0.0624. Blum et al. introduced an 
additional term to take boundary effects into account, 



(6) 



where n is the quantum number, which in our case is 
given by the Weyl number n-yveyi (see Eq. ([3])). In Fig. [3] 
number of nodal domains is shown versus the Weyl num- 
ber for the real (crosses) and imaginary (diamonds) part 
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FIG. 4: (color online) The variance of the number of nodal FIG. 5: (color online) The distribution of normalized nodal 
domains u n versus the Weyl number nwoyi for the real (blue domain areas Pa(s / s m in) is shown for the real (solid his- 
boxes) and imaginary part (red crosses). The blue dashed line togram) and the imaginary part (red dotted histogram). Ad- 
corresponds to the theoretical prediction. ditionally the expected algebraic decay with an exponent of 

r = 187/91 is shown as blue dashed line. 



of the wavefunction. The dashed line corresponds to the 
slope given by Ref. 1. The data were fitted with Eq. © 
and we found a=0.0594 and 6=1.231 for the real and 
a=0.0599 and 6=1.300 for the imaginary part. The value 
for a is in accordance with the expected slope of 0.0624. 

The boundary correction parameter 6 should scale with 
the ratio S/ vl, where S and A are the circumference and 
area of the billiard, respectively. The rescaled quantity 
b c = b^/A/S should be universal. In the present case 
we find values of 0.29 and 0.31 for 6 C for the real and 
the imaginary part, respectively. For two other billiards 
investigated experimentally,namely the half rough circle 
billiard with Q and without [|[ teflon insert, values of 
0.27, and 0.20 have been obtained for 6 C . The 6 C obtained 
for the billiard with a half circle teflon insert have to 
be taken with care, as the internal boundary between 
teflon and air and the integrable shape of the teflon insert 
may give rise to deviations. In case of the numerical 
calculations presented in Ref. d values of 0.25 and 0.21 
have been obtained for the Sinai and stadium billiard, 
respectively. In all cases the estimated error is of about 
10%. Thus b c is comparable for all billiards. 

For the variance a 2 of the number of nodal domains 
the percolation model predicts a 2 w 0.05rt 4]. In Fig. [4] 
the scaled variance cr 2 / n Weyi is plotted versus the Weyl 
number nwoyi for the real (blue boxes) and imaginary 
(red crosses) part. The variances have been calculated 
from the corresponding fits mentioned before and the av- 
erage has been performed over 30 consecutive wavefunc- 
tions. A good agreement is found. The larger variations 
at small Weyl numbers for the imaginary part are due to 
experimental noise: As in this regime the phase rigidity is 
often small, i. e., the imaginary part is small, noise in the 
experimental data can create small nodal domains par- 
ticularly close to the billiard boundary, leading to higher 
fluctuations. 

Another quantity of interest is the distribution Pa of 
nodal domain areas s. The area distribution has also 



been calculated in Ref. H and is given asymptotically by 

P A (s) oc (—) (7) 

where r = 187/91 and s m i n = n(ji/k) 2 , the smallest pos- 
sible area for fixed wavenumber k, where ji is the first 
zero of the Bessel function Jq(x). In Fig. [5] the distribu- 
tion of normalized nodal domain areas is shown for the 
real and imaginary part. All evaluated wavefunctions 
have been used to obtain the distribution. Additionally 
the expected slope of the algebraic decay with r = 187/91 
is shown as a dashed line and is in agreement with the 
experimental histograms. Deviations at large values of 
■s/smin > 80 are due to poor statistics. Here only a few 
entries per bin are available. 

No difference between the expected behavior for closed 
chaotic systems and the behavior for the real and imag- 
inary part of the wavefunction in open chaotic systems 
is found. A global phase rotation had been performed to 
uncorrelate real and imaginary parts, but the same be- 
havior is found for the number of nodal domains (Fig. [3]), 
its variance (Fig. [4| and the area distribution (Fig. [5]), 
using the original data without any global phase rotation. 

III. NODAL DOMAIN DEPENDENCE ON 
GLOBAL PHASE 

In the preceding section we have mostly neglected the 
fact that the billiard is open and treated the real and 
imaginary part just like a real valued wavefunction in 
case of closed systems. But as the billiard is open we 
have additional parameters like the rigidity \p 2 \ or global 
phases ip g or even new quantities like currents j or vor- 
ticities u> [lflU^]- In Fig. [6] we present the nodal domains 
and the number of nodal domains of the real part of a 
wavefunction at Weyl number n-\y ey i ~ 223 with a phase 
rigidity \p\ 2 0.81. At the top the number of nodal 




FIG. 6: (color online) (Upper part) Number of nodal domains i/ n as a function of a global phase <p g for a wavefunction at 
^Wcyi ~ 223. Marked are the phases <p g , for which the nodal domains are shown below, corresponding to <p> g /n — 0, 0.225, 
0.35, 0.475, 0.5, and 0.5525. Areas, where appearances and disappearances of nodal domains can be seen, are highlighted. The 
inserts show the corresponding ipi versus ipn plots, see also Fig. [2] 



domains for the real part is shown as a function of the 
global phase ip g . For ip g — real and imaginary part 
are uncorrelated and for tp g = n/2 real and imaginary 
part have changed their identities and are uncorrelated 
again. In the lower part a selection of nodal domains is 
presented for phases indicated by arrows in the upper fig- 
ure. Additionally the ipj versus tpn plots are shown, see 
Fig. [5J Highlighted rectangles emphasize regions, where 
changes of the number of nodal domains occur, due to re- 
arrangements of nodal lines. While the phase is changing 
the nodal lines are shifted and permanently dissolved and 
reconnected. The only points fixed are the nodal points 



of the corresponding modulus |V>| 2 , corresponding to el- 
liptic points in the flow. Nodal lines can only intersect at 
the saddle points corresponding to hyperbolic points in 
the flow. Their position is also independent of the global 
phase (see Ref. [29| for a more detailed discussion of these 
features) . 

Because the phase rigidity for this wavefunction is 
quite large (|/9| 2 =0.81) the modulus of the imaginary 
part \tpi\ is much smaller than the modulus of the real 
part Therefore a small change in the global phase 

(p g i ► (p g ~f" dip for (p g close to zero (or close to 7r) has only 
a weak effect on the real part of the wavefunction and on 
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FIG. 7: Number of nodal domains v n as a function of a global 
phase ip g for a wavefunction at n-wcyi ~ 229. 




FIG. 8: (color online) Correlation of the signed areas as a 
function of the global phase <p g for three different wavefunc- 
tions at Weyl numbers nwcyi ~ 229, 222, and 223 with rigidi- 
ties /9=0.04, 0.48, and 0.81, respectively. The curve with the 
smallest curvature belong to p =0.04, and the one with the 
largest curvature belong to p =0.81. The dotted lines show 
the results of Eq. (|12|l • The first and the last correlation cor- 
responds to the wavefunctions used to produce Figs. [6] and 

m 



its nodal pattern. Here the number of nodal domains 
v n is less likely to change as a function of phase. On 
the other side for ip g close to tt/2, there is an increased 
probability for changes in the nodal pattern. Indeed, if 
the global phase approaches 7r/2, changes of v n become 
more frequent. In Fig. [7] the number of nodal domains 
as a function of the global phase is shown for wavefunc- 
tions at Weyl numbers n\v e yi=229 with |p| 2 =0.04. In this 
case the phase rigidity is small, i. e. real and imaginary 
parts are of the same order. Therefore already for small 
global phases the number of nodal domains changes and 
a suppression of changes in v n for small phases is not ob- 
served - instead the changes occur quite uniformly over 
the complete range of global phases. This already shows 
that the phase rigidity can be observed from the number 
of nodal domains, as a function of the global phase, at 
least qualitatively. 

It is a non-trivial problem to obtain an expression for 



the number of nodal domains as a function of the global 
phase based on the random-wave (or the percolation) 
model. However there are other quantities which describe 
the nodal pattern as a function of the global phase and 
which are less resistent to a theoretical approach. Any of 
these is suited to exhibit the interplay between openness 
and nodal domains. One possibility is the global phase 
autocorrelation of the signed areas of the wavefunctions 
defined by 

C(tf) = (sgnCV'R.^i) ' S gn(-0R, V 2)>; <P = <P2-<Pi, (8) 

where sgn(x) denotes the sign of x, and (•) denote an 
average over the area. ipn. v denotes the real part of the 
wavefunction for the value <p of the global phase, ipi is 
chosen to be the global phase, where the real and imagi- 
nary parts are uncorrelated, in our case p\ = 0. Equiv- 
alently one can write 



C(<p) 



A, 
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where A is the area of the billiard, A$ t+ is the area where 
V'fl.^i an d ipR,<pi have the same sign, and A#,_ = A — 
A$ t+ is the area where they have opposite signs. Thus 
the signed area correlator measures the fraction of the 
billiard area which changes sign as i? increases. 

We can give an expression for the autocorrelation func- 
tion C(ip) based on the random wave model. At a given 
point we may write the wavefunction at tp = as 



V^o = if>R + itpj =xi + i£x 2 



(10) 



where x\ and x 2 are independent and equally distributed 
Gaussian random variables, and £ is a real constant, that 
is related to the rigidity via 



m -mi 



(11) 



In general, random wave models of the given type only 
define the rigidity on the mean over all realizations (nei- 
ther the norm of the real part nor the norm of the imag- 
inary part are fixed). For the present the fluctuations 
in ij) 2 R and ipj just reflect the fact that we describe one 
point in a wavefunction while the rigidity is obtained by 
an integral (space average) over the billiard. Other ap- 
plications of random wave models may need to take more 
care and define a model where the rigidity does not fluc- 
tuate from one realization to another. 

Within the above variant of the random wave model 
the autocorrelator can easily be calculated as 



C(<p) 



1 



dx\dx2 e Xl 



xsgn(:ri)sgn(:Ei cos(<^) + £x 2 sin(y>)) 
f — sgn( cos a)sgn( cos a cos ip + £ sin a sin tp) 

JO 27T 

(12) 
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FIG. 9: (color online) Left the probability current density j 
is shown for the whole billiard at 13.84 GHz. On the right 
side the corresponding vorticity uj is plotted. The lower part 
of the figure shows a magnification of the region marked by 
squares in the upper figure. In the zoomed figures vortices and 
saddles are marked by circles and crosses. Open circle denote 
clockwise and filled circles denote counterclockwise vortices. 

In Fig. [5] the global phase autocorrelation of the signed 
areas is shown for three different wavefunctions at ap- 
proximately the same Weyl number nwcyi ~ 229, 222, 
and 223 but for very different phase rigidities |p| 2 =0.04, 
0.48, and 0.810, respectively. Additionally, the results 
from Eq. (fT2"]) are plotted as dotted lines. An excellent 
agreement between experiment and theory is found, es- 
pecially as there is no free parameter, since the only pa- 
rameter, the phase rigidity, has been determined directly 
from the wavefunctions. 



IV. NODAL DOMAINS FOR THE VORTICITY 

Since the billiard is open there is a flow through the 
system which, in the electromagnetic case, is described by 
the Poynting vector [301 ] - Quantum mechanically it cor- 
responds to the probability current density and is given 
by 

j oc Im[V>*(r)VV>(r)] . (13) 

In Fig. [9] (left) the probability current density j is 
shown for the same wavefunction as used in Fig. [SJ A 
complex flow structure is observed with numerous ellip- 
tic fixed points (vortices) and hyperbolic points (sad- 
dles). The vortices and saddles are marked by circles 
and crosses in the lower part of the figure. 

Another quantity of interest is the vorticity given by 

W = (V^flXVj,^/) - (V„Vfl)(V*^/). (14) 
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FIG. 10: (color online) The number of nodal domains v n ver- 
sus the Weyl number nwcyi for the vorticity to (diamonds). 
The red dashed line correspond to the theoretical prediction 
v n = 0.0624n of the percolation model UJ. The blue solid 
lines are fits including boundary effects with Eq. with 
a linear slope a = 0.161. 



which, up to a constant factor, is just the curl of flow. It is 
shown Fig. [5] on the right hand side. For more details on 
distributions and correlations of the probability current 
or vorticity we refer to Refs. [TH, [28|. One observes clearly 
a nodal line pattern for the vorticity. The nodal lines of 
the vorticity correspond to the unstable manifolds of the 
dynamics. Accordingly, nodal lines of the vorticity will 
always cross exactly at the saddle points and each vortex 
has a separated nodal domain with an area depending 
on the vortex strength and the distance to other vortices 
and saddles. A type of irregular checkerboard pattern of 
two saddles and two vortices, one clockwise, the other one 
counterclockwise, with a nodal line crossing between the 
four of them is created and can be seen in Fig. [9](upper 
right). Due to measurement errors but also already due 
to the discretisation and the bilinear interpolation these 
crossings will be turned into avoided crossings, thus con- 
necting the nodal domains more or less randomly. This 
corresponds to the spirit of percolation theory. 

In Fig.[Tn]we present the number of nodal domains for 
the vorticity. A fit with Eq. © yields a linear slope of 
a=0.161 and 6=1.65. Therefore we find again a depen- 
dency like in Eq. ©. The slope is larger than the slope 
expected in the case of real wavefunctions (0.0602). Pre- 
dictions for the linear increase of the vorticity in terms 
of a percolation model are still missing. 



V. SUMMARY 

To summarize, we have shown in this paper that the 
nodal domains of the real and imaginary part of the wave- 
functions in open systems behave like the nodal domains 
for the wavefunction in closed systems. Additionally we 
have calculated the global phase autocorrelation of the 
signed areas as a function of the global phase. This quan- 
tity was shown to be an indicator of openness, i.e., the 
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phase rigidity. An effect of the rigidity is also present in 
the fluctuations of the number of nodal domains as the 
global phase is varied, though an analytical description 
is beyond present knowledge. Nodal domains in open 
systems can also be defined for other quantities like the 
vorticity, which we have shown here. They also seem to 
have a linear behavior (plus square root corrections) but 
with a different slope in comparison to the number of 
nodal domains in wavefunctions in closed systems. 
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